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Abstract. As a quantum device, a quantum heat engine (QHE) is described by a 
Hermitian Hamiltonian. However, since it is an open system, reservoirs have to be imposed 
phenomenologically without any description in the context of quantum mechanics. A non- 
Hermitian system is expected to describe an open system which exchanges energy and particles 
with external reservoirs. Correspondingly, such an exchange can be adiabatic in the context 
of quantum mechanics. We first propose a non-Hermitian QHE by a concrete simple two-level 
system, which is an S' = 1/2 spin in a complex external magnetic field. The non-Hermitian 'P'T- 
symmetric Hamiltonian, as a self-contained one, describes both working medium and reservoirs. 
A heat-engine cycle is composed of completely quantum adiabatic processes. Surprisingly, the 
heat efficiency is obtained to be the same as that of Hermitian quantum Otto cycle. A classical 
analogue of this scheme is also presented. Our finding paves the way for revealing what the role 
of a non-Hermitian Hamiltonian plays in physics. 
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1. Introduction 

Non-Hermitian system has received intensive studies due to two reasons: It can possess full real 
eigenvalues [1] which is fundamental in the quantum world. A non-Hermitian potential breaks 
the probability conservation, which is supposed to characterize the source or drain of particles and 
energy. The connection between the imaginary potentials and the environment has been investigated 
in discrete system [5] . A non-Hermitian system could be employed to describe an open system which 
exchanges energy and particles with external reservoirs. An alternative way to deeply understand 
the exact meaning of a non-Hermitian Hamiltonian in physics is to find out its function which is 
similar to some phenomena in practice. An exemplified example is the invention that an imaginary 
potential can be the function of laser or anti-laser [3l|4] when spectral singularity is reached [5]. 

The topic of the quantum heat engine (QHE) [6l [71 |8l 0 [101 El lEl US] has attracted a 
lot of interest since it was first proposed by Scovil and Schultz-Dubois M- A QHE converts 
heat into useful mechanical work from quantum working mediums, such as multilevel systems m. 
harmonic oscillator systems [ISllIillIlllIl], spins or coupled spins |2nj , optomechanical systems |21] , 
relativistic particles [22j and so on. Many investigations have been carried out to explore various 
possible improvement of the heat engine that the efficiency of a QHE can surpass the efhciency of a 
classical Carnot heat engine, and a QHE can extract work from a single heat bath . Due to the 
quantum properties of the working medium, or the effects of the quantum heat bath, some unusual 
and exotic phenomena have been manifested by considering and using a squeezed reservoir [241125] , 
quantum coherence [23], or coupled spins [26l[271 [JHIES] . Moreover, via constructing a QHE which 
is a two-level quantum system and undergoes quantum adiabatic process and energy exchanges 
with heat baths at different stages in a work cycle, some important aspects of the second law of 
thermodynamics have been clarified by Kieu [30] . Very recently, many other investigations have been 
carried out and demonstrated about the Carnot statement of the second law of thermodynamics and 
the quantum Jarzynski equality in quantum systems described by pseudo-Hermitian Hamiltonians 

m- 

While, in contrast to those quantum devices, external reservoirs are inevitable. QHEs 
considered in the literature are mainly described by a Hermitian Hamiltonian and imposed 
reservoirs. With the proposal and development of theoretical explorations [32l [33l [34], non- 
Hermitian quantum theories arise as effective descriptions of certain open quantum systems, and 
therefore induce to an effective non-Hermitian QHE in the presence of absorption and gain. It is 
natural to establish a non-Hermitian description for QHEs. 

In this paper, we employ a non-Hermitian PT-symmetric Hamiltonian to study a QHE. It is 
a two-level system, which describes an S = 1/2 spin in a complex external magnetic field. The 
non-Hermitian QHE behaves quite differently from a Hermitian one in all the approaches so far, 
providing an alternative description. We do not provide a Hamiltonian to describe the working 
medium, whereas a fully description of working medium and external reservoirs is included. The 
process of exchanging heat can be adiabatic in the context of quantum mechanics. We will show that 
such a non-Hermitian system can fully describe a QHE without imposed external reservoirs, which 
can operate at Otto efficiency for an optimal cycle. In order to get a clear physical picture of the 
non-Hermitian engine we construct a classical analogue of this scheme, which is also a variable-mass 
cycle. 
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Figure 1. (Color online) Schematic diagrams of a variable-mass Otto cycle for the non-Hermitian 
QHE. The cycle is shown as a closed loop in the (a) J — (f> and (b) e+ — P+ planes, respectively, 
consisting of four steps, two isospectrum processes A B and C D; and two adiabatic 

processes B C and D —> A. According to Eq. II17II . the mechanical work done by the engine 
equals the surface integral of a nonzero function over the surface enclosed by the loop. 


2. Model 


We consider a non-Hermitian spin-1/2 system in an external magnetic field, which can be described 
by the following Hamiltonian 

H = B-a, (1) 

where B = (Bx, By, i^J ,) is a time-dependent complex magnetic field and 3 = {(Jx, (TyjO-z) is Pauli 
matrices. Taking 

Bx = J cos (j), By =—J sin (j), (2) 

cos 6 = i'-^l\/\ — 7 ^, 


we rewrite the Hamiltonian as the form 


H = J{t) v'l-y" 


cos 9 

sinOe-^^'dl 


sin6»e*‘^(‘) \ 

— cos 9 J 


(3) 


where J (t) and are real functions of time, and 9{'y) is a complex (real) constant. Obviously, 
the Hamiltonian H is PT-symmetric, i.e., [P,H] ^ 0 and \T,H] ^ 0, but \PT,B\ = 0, in the 
sense of BiT = —i and V = Ox- The diagonal elements of the matrix can be regarded as imaginary 
on-site potentials of a two-site tight-binding model, the existence of which violates the law of 
conservation of mass and energy. It is thought as an open system with the source and drain. The 
aim of this paper is trying to establish a full quantum mechanical description of the QHE. We will 
consider a non-Hermitian two-level quantum engine, in which the source and drain are thought as 
the channel to the heat bath. To this end, we will seek the solution of the Schrodinger equation 
of the system. Considering the adiabatic time-evolution of an initial eigenstate of B (0) under the 
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quantum adiabatic condition, we have 
I'hA it)) = it)), 

where \tf}\ (<)) (A = ±) is the instantaneous eigenstate of H {t), 


1 -^+) = 


cos I 
sin 


>IV'-) = 


— sin ^ 
cosfe-*'^ 


(4) 

(5) 


with eigenvalue e\ = XJ (t) v^l — 7 ^ and J (t) > 0 (see Appendix). Phase Ax is the key quantity 
of this work, which can be expressed as 

Aa = - Xy/l - 7 ^ [ J (t) dt 
Jo 

+ ^[^W-^(0)] (l-zAy/v'l-y^) . (6) 

In the rest of the paper, our investigation does not involve a coherent superposition of \^px (t)). We 
neglect the real part of Aa and take 

AA = -^[<^(t)-<(>(0)]/v'W. (7) 

We note that Aa is path independent and proportional to the difference of </). And the nonzero 
imaginary adiabatic phase is quite essential to the thermal process of QHE, since it can vary the 
amplitude of instantaneous eigenstates [35) . leading to the change of population distributions. The 
non-Hermitian Hamiltonian is rather different from a Hermitian one, treating the system as a whole 
in describing the working medium and reservoirs. 

Before start dealing with a heat engine cycle, we would like to briefly remark on the time 
evolution of a non-Hermitian system, which involves both non-unitary evolution and complex 
geometrical phase of a state. The Schrodinger equation is the common basis of both conventional 
and non-Hermitian systems. A quantum state evolves along the solution of the Schrodinger 
equation. According to the solution, the Dirac probability is no longer conservative. In this sense, 
the non-unitary evolution and the complex geometric phase are nothing but the direct and natural 
results of the Schrodinger equation. 


3. All-quantum-adiabatic cycle 

Now we investigate a thermodynamic process for an adiabatic passage of the non-Hermitian system. 
For a traditional (Hermitian) QHE, the temperature of a quantum working medium is introduced 
based on the thermalization assumption [SlIlolEolEe], that is, the working medium has the same 
distribution of occupation probability as the heat bath in thermodynamic equilibrium. Then for 
such a small quantum system, the probability distribution is related to a parameter called spectral 
temperature. In this paper, we still employ such a concept for the non-Hermitian two-level system, 
as an assumption. Consider an initially prepared mixed state, which has Dirac probability P\{t = 0) 
on the eigenstate I^/’a {t = 0)) of H (0) with eigenenergy AJo\/l — 7 ^. The system parameters are 
taken as (Jo,(()o) initially where Jq = J (0) and (po = (p^O). And the initial populations for the 
upper and lower levels are T+(0) = Po and T’-(O) = 1 — po, respectively. We express this mixed 
state by a density matrix 

p(0) = ^ Pa(0) I^A (0)) (^A (0)1, 

A=± 


( 8 ) 
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Table 1. State parameters for a quantum variable-mass Otto cycle illustrated in Fig[T] Here 
T, e+ and U are in the unit of — 7 ^. 



J 

0 

T 

£+ 

P+,P- 

u 

s 

A 

Jl 


2 Jl 

Jl 

PoA-Po 

Jl (2po - 1) 

-k-BPo Inpo 

fcs ln(pg'^-l) 

-fee (1 -po)ln(l -po) 

B 


4>2 

2 Jl 

Jl 

PoA (1 

po (e + r')^i 

-feppoCln (PoC) 

fcs ln[(po 

-fee (1 - Po) In [(1 - Po) 

C 

J2 

4>2 

2 J 2 

J2 

PoA (1 

Po (C + J2 

-r'^2 

-^BPo^ln {poO 

fcs ln[(Po 

-fee (1 - Po) In [(1 - Po) 

D 

J2 

</>l 

2 J 2 

J2 

PoA-Po 

J2 (2po - 1) 

-fcePo Inpo 

fcB ln(pp “^-l) 

-fce (1 -po)ln(l -po) 


Table 2. Variances of state parameters related to four processes in a quan tum variable-mass 
Otto cycle illustrated in Fig[T] Here AT and AU are in the unit of \/l — 7 ^. 



AT 

AU 

AQ 

AW 


AS 


A^B 

4 Jl In C 

Po (C + - 2 ) Jl 

+ (1 - r') Jl 

AU 

0 

-/cBPolnfpo - fce (1-Po) 

fen ln[(p-i-l)C-2j In(p-^-l) 

X In 



B^C 

2(J2-Ji) 

Po iJ 2 - Jl) 

0 

AU 


0 


fen ln[(p“^-l)4-2j 

-r' iJ2 - Jl) 



C^D 

4J2lnf 

-po (C + - 2) J 2 

-(i-rJJ2 

AU 

0 

fcaPo In (po + A:b (1 - Po) 

ks ln[(p“'^-l)4-2j ln(p"'^-l) 

X In 



D ^ A 

2(Ji-J2) 

(Jl - J 2 ) (2po - 1) 

0 

AU 


0 


ks ln(pg'^-l) 




where eigenstate \ijj\ (0)) is normalized in the context of Dirac inner product (see Appendix). Here 
we assume that the Boltzmann distribution still refers to Dirac probability even for a non-Hermitian 
system. When a two-level system couples to a heat bath, both systems obey the thermal equilibrium 
Boltzmann distributions and then the two-level mixed state has its temperature. A density matrix 
would represent a thermal equilibrium state with temperature 

^ 2Jov^l - 7^ 

° fee In (pq ^ - 1 ) ’ 


(9) 
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where is the Boltzmann constant. In quantum statistics, the definition of temperature does not 
require the preservation of particle probability. It is only determined by the distribution of particle 
probability on each level. For a two-level system, one can always find a temperature which only 
depends on the ratio of the probabilities and the energy difference of two levels, i.e., 


^ _ e+ - e_ 
P_ knT 


( 10 ) 


where P\ and £\ are the particle probability and the energy of two levels. At time t, state |'0 a (0)) 
evolves into state jd/A (t)) under the time-dependent Hamiltonian H (t). We note that the initial 
state does not involve a coherent superposition of IV'-i- (0)) and IV'- (0)). Then the evolved mixed 
state does not involve a coherent superposition of (t)) and Idf- (t)). We only concern about 
the time evolution of the pure state IdtA (t)). 

It is an acceptable consensus that the Schrodinger equation is the common basis for both 
Hermitian and non-Hermitian quantum mechanics. Whether the operator H (t) is Hermitian or 
not, an evolved state, vector |4 'a (t)) should be the solution of the Schrodinger equation 

*||vhA(0)=^WI^A(t)), (11) 

with the initial condition |'I'a (0)) = ji/'A (0)). Accordingly, the density matrix of the evolved mixed 
state becomes 


p(t) = ^PA(0)|d/A(t))(d'A (t)|, 
A=± 


( 12 ) 


where |'I'a (t)) is probably not Dirac normalized so as to contribute to the population Px{t). 
However, it is tough to get an analytical solution for an arbitrary time-dependent operator H (t). 
Fortunately, when we only consider the time evolution under the quantum adiabatic condition, 
i.e., there is no transition between two levels, an approximate solution can be obtained. Then the 
density matrix of the evolved mixed state can be expressed as 


P{t) = 


A=± 


where 


P\{t) = Pa(0) exp[-2Im(AA)]. 

And the corresponding temperature of the thermal state p(t) is 

2Jit)x/T^ 


T(t) = 


kB In [(Po^ - 1 )^"^] ’ 


(13) 


(14) 


(15) 


where ^ = expjy {(j>{t) — /-^/l — 7 ^} characterizes the amplification or attenuation of the whole 

probability in the working medium. 

Under the quantum adiabatic condition, the solution of Eq. dU tells us that any evolution along 
an arbitrary path in the J — (p plane is always adiabatic in the context of quantum mechanics, i.e., 
there is no transition between two levels. Nevertheless, it does not mean that the occupation 
probabilities for the levels and their distributions are always remaining invariant due to the 
contribution of the Aa in Eq. 0 , which are unusual and exotic properties in the non-Hermitian 
system in contrast to the Hermitian one. Moreover, we would demonstrate that the “adiabatic” in 
our paper for the non-Hermitian Hamiltonian may be different from the one in the Hermitian system 
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because the occupation probability for each level is not invariant due to the image part of A^. But 
what both adiabatic processes of the Hermitian and non-Hermitian system have in common is that 
there is no tunneling between each level. Therefore here we could call it “generalized adiabatic” 
in our non-Hermitian system. Since all the processes considered in our work are “generalized 
adiabatic”, we only write it as the “adiabatic” for simplification. Now we consider a loop in the 
J — (j) plane as a thermal cycle. 

According to the first law of thermodynamics, the total heat transferred in the cycle is 

AQ = i dQ = (1)'^ e\dP\ 

Jc Jc ^ 


27 




1 


= — // cosh In W-1 —lnMdJd(/), 


Po 


(16) 


where Z = 2cosh(ln^^ — 1) is the partition function. Here A denotes the area enclosed by the 
loop c in the J — (j) plane. We can see that AQ is nonzero for a non-trivial cycle, since a hyperbolic 
cosine function is positive definite. Similarly, the net work done in the cycle is 

AW = /dW = / ^ PAdsA =-AQ, (17) 

Jc Jc ^ 


which leads to 


AU = AQ + AW = 0. (18) 

Furthermore, a straightforward derivation shows that the entropy variation in the loop is 

^ fc ^ i- ^ ^ ^ 

According with the traditional thermodynamics, these indicate that the internal energy U and 
entropy S are state functions, which are only determined by the system parameters {J,(l)). We 
would like to point out that the facts of AQ ^ 0 and AW ^ 0 in the loop are exclusive for a 
non-Hermitian system. Taking 7 = 0 , the system is reduced to a Hermitian one. We can see that 
AQ = AW = 0, which indicates that the Hermitian system is trivial as a QHE. More detailed, 
as mentioned above, the complex parameter represents the effect of reservoirs phenomenologically. 
Once 7 = 0 , the working medium is decoupled from reservoirs. There is no exchange between the 
working medium and the reservoirs. Then the heat engine becomes trivial, since we do not assume 
extra reservoirs imposed on the system. This is different from the conventional Hermitian QHE, 
which is always accompanied by the reservoirs. This reflects the key point of our work, replacing the 
imposed reservoirs by the imaginary parameter 7 . In the context of quantum mechanics, processes 
for an arbitrary loop is spontaneously adiabatic without any extra condition. This may enhance 
the feasibility of experimental realizations of the present scheme. 


4. Cycle with Otto efficiency 

So far we have established an alternative description for QHEs via a non-Hermitian Hamiltonian. 
Remarkably, all the thermal processes are adiabatic in a quantum mechanics manner. By this we 
mean that irrespective of whether there is energy transfer between the heat engine and reservoirs 
or not all the processes are adiabatic evolutions driven by a time-dependent Hamiltonian. To 
understand such quantum processes, it is useful to contrast them to classical thermal processes. 
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Figure 2. (Color online) A classical analogue to the quantum variable-mass cycle consisting of 
four steps, two isothermal processes and two adiabatic processes. Here the isothermal process is 
associated with particle transfer. The cycle is shown as a closed loop in the (a) P — V and (b) 
T — M planes, respectively. 


In the following, we try to connect such a description to a classical engine, showing that the non- 
Hermitian QHE is not a non-physics toy model. 

To this end, we consider a specific cycle, which consists of four processes connecting four points, 
A: B: {Ji,(j) 2 ),C: (J 2 , and £>: (J 2 ,^!)i). There are two types of processes involved: (I) 

Isospectrum process, A ^ B and C ^ D\ (II) Adiabatic process, B ^ C and D ^ A. A schematic 
illustration is listed in Figs. [iDa) and[T](b). In type I process, varying (p can change the populations 
on two levels due to the imaginary adiabatic phase Ay, but leave the level structure unchanged. The 
level structure can be completely characterized by quantity £+. While, in type II process, varying 
J can change the level structure e+, whereas leave populations on two levels unchanged. Based on 
the adiabatic solution of Eq. o, quantities, including the state parameters at each point and their 
variances in each process, are obtained and listed in Tables [T] and [U respectively. According to the 
definition in the Refs. [3011II0], our calculations indicate that such an engine can operate at Otto 
efficiency 


77 = 


Qi + Q 2 

Q2 




( 20 ) 


Such a non-Hermitian QHE does not violate the Carnot efficiency limit which is of great 
significance in thermodynamics. This indicates that the non-Hermitian Hamiltonian still obeys 
some underlying rules in nature. 

We note that the probability is not conservative in the cycle, which differs from that in a 
Hermitian system. We refer to this cycle as a variable-mass Otto cycle. However, we would like 
to point out that the particle-number conservation still holds in such a cycle m- Here we use the 
term variable-mass in order to correspond with its classical analogue. For a particle-exchange heat 
engine, it has been investigated in Ref. |38] . Moreover, to understand such a non-Hermitian cycle 
and its related issues, we propose a variable-mass cycle for classical ideal gas in the following. 
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5. Classical analogue 

In order to obtain a clear physical picture of the non-Hermitian QHE we construct a classical 
analogue of this scheme, which is a variable-mass cycle and consists of the following four steps, 
which is schematically illustrated in Fig. [2] 

1. Reversible isothermal addition of gas at the high temperature Ti. During this step a volume 
AVi of the gas with the same temperature Ti and the same density as the working medium is 
added. It does no work on the surroundings. The temperature of the gas does not change during 
the process, and thus the addition is isothermal. The gas addition is associated by absorption of 
heat energy Qi and increase of entropy AS” = QifTi from the high temperature reservoir. 

2. Reversible adiabatic expansion of gas. For this step, the gas is allowed to expand and does 
work on the surroundings, whose expansion causes it to cool to the low temperature T 2 , with the 
entropy remaining unchanged. 

3. Reversible isothermal deduction of gas at the low temperature T 2 . During this step a 
volume AV 2 of the gas with the same temperature T 2 and the same density as the working medium 
is removed. Such an amount of the gas contains heat energy Q 2 and entropy AS” = Q 2 /T 2 , which 
flows to the low temperature reservoir. 

4. Reversible adiabatic compression of gas. During this step, the surroundings do work on the 
gas, increasing its internal energy and compressing it, causing the temperature to rise to Ti. The 
gas goes back to the state at the start. 

For a standard cycle, mass is conserved, i.e., the working medium does not lose or increase mass. 
The present cycle transfers mass from hot to cool reservoirs, associating with the heat transfer. The 
essence of this cycle is that there is a volume of gas transferred from hot to cool reservoirs, rather 
than transfer the heat only. Then the heat efficiency is always 1 — T 2 /T 1 , which is independent of 
the mass transferred. Therefore, the proposed variable-mass cycle would help us to make a better 
understanding of the non-Hermitian cycle in some sense. 

6. Summary and discussion 

We have studied a time-dependent non-Hermitian PT-symmetric two-level system, which could be 
regarded as a QHE. A non-Hermitian Hamiltonian turns out to be quite different from a Hermitian 
one. It does not solely describe the working medium, but offers a full description including the 
working medium and external reservoirs. We have shown that any cycle can be implemented by an 
adiabatic time-evolution along a quantum adiabatic passage. As an example, a specific cycle which 
operates at Otto efficiency is presented. Moreover, we have also discussed a classical analogue of 
this scheme, which corresponds to a variable-mass cycle, transferring not only heat but also mass 
from hot to cool reservoirs. Our main conclusion is that the proposed non-Hermitian Hamiltonian 
can provide an alternative description for the QHE, leading to a deeper and better understanding 
of the role that a non-Hermitian Hamiltonian takes in physics. 

Finally, we would like to point that there is a long way to go for connecting the current non- 
Hermitian 2-level system to a Hermitian description of QHE by using a rigorous mathematical 
formulation. It is due to the subtle relation between Hermitian and non-Hermitian Hamiltonians: 
on the one hand, a non-Hermitian Hamiltonian is usually thought to be a better candidate for 
describing some natural processes in an open system, which do not follow the conservation of law 
of mass and energy. However, on the other hand, it is still a challenge to get a perfect connection 
between a non-Hermitian Hamiltonian and a Hermitian one with a clear physical picture, although 
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7. Appendix 


The instantaneous eigenstates of H {t) and H'^ {t) with eigenvalue £a = A J (t) are 



cos I 
sin 16“"^ 


) .IV'-) = 



( 21 ) 


and 



( 22 ) 


respectively, which can construct biorthonormal complete set {IV'a) ; |?7a)} (A = ±), i.e., 


ivy I^a) = <5aa',X! 


(23) 


A 


Here we demonstrate that when computing the density matrix or detecting the probability in 
experiments, both eigenstates of H (t) must be normalized in Dirac inner product, i.e., we need to 
take |'(/'a) —>■ IV’a) /a/© with 0 = |cos|| + |sin|| . Now we extend the adiabatic theorem to the 
non-Hermitian Hamiltonian and we will show that such an extension is correct. This extension can 
be simply applied by taking the biorthonormal inner product to replace the Dirac inner product. 
Then under the quantum adiabatic condition 





{s+-e-f 


4J(t)(l-72) 


phase A\ can be expressed as 



(25) 


Similar works have been done in many systems [321 SHI SIl SH @3] , and this conclusion can also be 
obtained by an exact solution of Schrddinger equation 


Id/) = Id/), 


(26) 


in this concrete example. Here our PT-symmetric Hamiltonian H could be used to describe an 
open system, which is different from [33] that is used to describe a closed system. Actually, the 
propagator of the Hamiltonian H 



( 27 ) 
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can be obtained as the form 
C/ii = 


U 22 = 


cos (Or) — i— sin (Or) 
cos (Or) + sin (Or) 






C/i2= -lsin(Or)e*-"/2^ 

C /21 = - ^ sin(Or) 

when we take (p (r) = cut. Here T is time-order operator and 


= ^‘j(t)dC,0 = iv'4 + A2, 


A = 217 -I- LO. 

The adiabatic condition in Eq. (l24ll reduces to 

< 1 , 


1-72 

under which the Taylor expansion gives 


(28) 


(29) 

(30) 

(31) 


— Or = 


- - 7^ 



ijuj 
1 — ^2 


r 



— 7 ^ J (t) dt 


i^UJT 

VT^' 


(32) 


This solution indicates that the adiabatic process exists under a certain condition, i.e., there is no 
tunnelling between instantaneous eigenstates during the time evolution. Then the corresponding 
geometrical phase can be obtained naturally. While, in the article |45j . dynamics near exceptional 
point is considered and it is natural that the adiabatic theorem cannot hold in that situation. The 
process cannot hold when the system closes to the exceptional point (or quantum phase transition 
point for breaking VT symmetry). Actually, similar behavior can occur in a Hermitian system, for 
example, an adiabatic process cannot be achieved near a quantum phase transition point. A simple 
demonstration is Landau-Zener formula. 

Moreover, in the biorthonormal formalism, compared with right ket |'k) which obeys the Eq. 
the left bra (<i)| follows the different evolution equation 

1$) = Ht 1 $), (33) 

and correspondingly the propagator 

U = T exp(—i f (t) dt) (34) 

Jo 


can be obtained as the form 

Uii = 


A* 


cos (HV) - *— sin (H*r) 


r/2 
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U22 = 


cos (0*r) + i— sin (IIV) 




Ui2= -^sin(r!V)e*-"/2^ 
U2i= -^sin(r!V)e-*“"/2^ 


then we can prove 

i/i'c/ = 1. 


(35) 

(36) 


As a result, it is nnitary according to the biorthonormal inner product. However, we would like 
to point out that the complex phase for a non-Hermitian system is a natural result (see references 
[Ml SQum sails]). The Berry phase and the unitarity of time evolution are compatible with each 
other. 

Next, explicit derivations of Eq. (IISl, (II3, and ra are given. Together with = po^, 
P_ = (1 — Po)^~^ (here ^ is only the function of (j)) and e\, we get the heat exchange, the work 
done, and the variation of entropy of any arbitrary path in the J — (j) plane as following. 


and 


AQ = j>dQ = exdPx 


- 7 ^ 


Jd 


Z 


Jd 


sinh (inW-1 — In^ 

V V Po 


= ‘^ j) J cosh ^In ^ — 1 — In^^ dcj) 

= ^ JJ cosh ^In ~ 1 ~ dJd^, 


AkE = j>dW = Pxdex 


dJ 


Z 




dJ 


= ■ T///“'■ 0" 


Obviously we have 

AQ + AIT = 0. 


(37) 


(38) 


( 39 ) 
































Non-Hermitian heat engine with all-quantum-adiabatic-process cycle 


12 


Similarly, the variation of entropy is 
AS = 


^ dg = (f 


T 


T 


A 




Z /e 

_ 2fcB7 

“ 

= 0 . 


Po 

In 


In M d 


sinh InW-1 — Inf 

Po 


Po 


— 1 — In ^ ) cosh I In 


Po 


-1-lnnd^ 


(40) 
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